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Mesoscopic models play an important role in our understanding of the deformation and flow 
of amorphous materials. One such description, based on the Shear Transformation Zone (STZ) 
theory, has recently been re-formulated within a non-equilibrium thermodynamics framework, and 
found to be consistent with it. We show here that a similar interpretation can be made for the 
Soft Glassy Rheology (SGR) model. Conceptually this means that the "noise temperature" x, 
proposed phenomenologically in the SGR model to control the dynamics of a set of slow mesoscopic 
degrees of freedom, can consistently be interpreted as their actual thermodynamic temperature. 
(Because such modes are slow to equilibrate, this generally does not coincide with the temperature 
of the fast degrees of freedom and/or heat bath.) If one chooses to make this interpretation, the 
thermodynamic framework significantly constrains extensions of the SGR approach to models in 
which a; is a dynamical variable. We assess in this light some such extensions recently proposed in 
the context of shear banding. 



I. INTRODUCTION 

Developing predictive theories for the deformation and 
flow of amorphous materials remains an ongoing chal- 
lenge. The class of materials concerned is large, includ- 
ing not only colloidal glasses, emulsions and foams but 
also molecular glass formers, metallic glasses and possi- 
bly granular materials. The development of such predic- 
tive theories would therefore be of great practical signifi- 
cance as well as fundamental theoretical interest ^] . Two 
key challenges in understanding deformation and flow in 
amorphous materials are the absence of a reference con- 
figuration relative to which strain can be defined (in con- 
trast to ordered crystalline solids), and the difficulty of 
encapsulating local flow events and/or larger-scale stress 
redistributions in a coarse-grained variable or continuum 
field (analogous to the dislocation density for crystals). 

First principles approaches to glass rheology have 
made significant progress but entail major approxi- 

mations (such as mode coupling theory), usually followed 
by further simplification 4]. In the absence of a com- 
prehensive microscopic theory, a number of models have 
been developed directly at a mesoscopic level to describe 
amorphous flow, see e.g. [5l-[l0|. Two such approaches 
that have been investigated in some detail are the shear 
transformation zone (STZ) approach, reviewed rece ntly 
in [ll|, and the soft glassy rheology (SGR) model fl^ - flSf . 
One important feature of the STZ theory is that it can 
be given an unambiguous thermodynamic interpretation, 
as discussed in detail in a series of recent papers [iGl - fTsI . 

Our purpose in this paper is to show that the SGR 
model can likewise be cast consistently within a non- 
equilibrium thermodynamic framework. This has two 



main benefits. Firstly, it gives an alternative way of 
looking at the effective temperature parameter x that 
was previously proposed within SGR to set the noise 
level for stochastic mesoscopic dynamics. This param- 
eter can now, if desired, by interpreted as a genuine 
non-equilibrium thermodynamic temperature, governing 
a subset of degrees of freedom whose dynamics causes the 
system to move among its various "inherent structures" 
or energy minima. Secondly, this interpretation signifi- 
cantly constrains how SGR can be extended to allow x to 
evolve in time and space: it suggests a specific form for 
the driving term in the dynamical equation for an evolv- 
ing X. We find that this driving term is consistent, in 
steady state, with one of two extensions of SGR postu- 
lated recently to understand shear banding in amorphous 
fiows [lH, but not with the other. Similar remarks ap- 
ply to two models for shear-thickening fluids developed 
previously [T9| . 
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We begin by reviewing briefly the SGR model, and 
then the non-equilibrium thermodynamic framework of 
Bouchbinder and Langer [T6l-[ll| (Sees. M and ini| . 
Sec. lIVI has the core of our argument, which shows that in 
this framework SGR can indeed be written as a thermo- 
dynamically consistent model, and that this consistency 
forces X to be the non-equilibrium thermodynamic tem- 
perature of a set of slow degrees of freedom. In Sec. |V] 
we discuss to what extent this argument can be made 
to cover also generalizations of the SGR model. Sec. IVII 
has a summary of our results and discusses the implica- 
tions for various extended SGR models that treat X 0(S Oil 
dynamical variable. 
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II. THE SOFT GLASSY RHEOLOGY MODEL 

The SGR model [ll-[l5l| describes a sample of amor- 
phous material as a collection of mesoscopic elements, 
chosen large enough that it makes sense to define local 
strain and stress variables. The local strain I is defined 
relative to a local equilibrium state 20], and is assumed 
to change in time by following the change in the macro- 
scopic shear strain 7, so that / = 7. We focus here on 
pure shear strain deformation but note that the model 
can also be extended to general deformation tensors [2l| . 
The local strain cannot grow indefinitely: once the stored 
elastic energy, which we write as lkve/2)P in terms of the 
local shear modulus k and the element volume We, be- 
comes close to a local yield energy E, a yield event will 
take place. This is assumed to reset I to zero and create a 
new local equilibrium configuration with correspondingly 
a new yield energy, drawn from a distribution p{E). The 
final ingredient is that this yield process is taken as ac- 
tivated by some effective temperature x. This is based 
on the intuition that every yield event elsewhere in the 
material causes a 'kick' locally via the associated stress 
redistribution, and that many such kicks add up to an 
effectively thermal noise. 

The resulting equation of motion for the distribution 
of strain and yield energy across elements, P{E, I, t), is 

P{E, I, t) = -^diP - T{E, l)P + Y{t)p{E, I) (I) 

Here V{E, I) is the yield rate from a given local config- 
uration, which represents the thermal activation in the 
form r(£;, /) = To exp[-(£; - kv^P/2)/x\. (We will leave 
T{E,l) general whenever possible in the following, to 
assess the scope of our arguments.) We use /cb = 1 
throughout, so that all temperatures are measured in 
energy units. The time-dependent total yield rate is 
Y{t) = J dEdlT{E,l)P{E,l,t). Finally, p{E,l) is the 
probability density of E and I for elements after a jump. 
This is equal to p{E)d{l) in the standard SGR model as 
described above. By choosing p{E) to have an exponen- 



tial tail 



-E/Eo 



one finds that a glass transition arises 
at x/Eq = 1. For x/Eq > 1 the material is ergodic (al- 
beit with infinite viscosity for 1 < x/Eq < 2) whereas for 
x/Eq < 1 the material is a nonergodic amorphous solid 
with time-dependent (aging) material properties 14]. 

The SGR shear stress is an average over elements. 



a{t) = J dEdlklP{E,l,t) 



(2) 



with k the local shear modulus as defined above. Start- 
ing from some initial condition P{E,l,0), this together 
with the evolution equation for P{E, l,t) in principle de- 
termines the stress a{t) for any imposed shear history 
7(t). 



III. NON-EQUILIBRIUM THERMODYNAMICS 
FRAMEWORK 



We next review the non-equilibrium thermodynamic 
framework of Bouchbinder and Langer (T6l - [T8j . The ba- 
sic premise is that the degrees of freedom of the ma- 
terial can be divided up, conceptually, into two weakly 
interacting subsystems. The degrees of freedom in the 
configurational subsystem (C) encode which of the expo- 
nentially many local energy minima, known in the glass 
community as inherent structures, the system finds itself 
in; this subsystem is 'slow' because at low temperature 
any motion between inherent structures has to be ther- 
mally activated or driven by external deformation. The 
kinetic-vibrational subsystem [K) gathers the remaining 
degrees of freedom that describe fast motion around these 
inherent structures. This subsystem is taken as strongly 
coupled to a thermal reservoir (i?) that sets the thermo- 
dynamic equilibrium temperature 9^. The total internal 
energy of the system is then written as 

C/tot = C/c(^c, A) + Uk{.Sk) + Ur{Sr) 

where Sc, Sk and Sr are the entropies of the config- 
urational and kinetic-vibrational subsystems and of the 
thermal reservoir. In Uci A denotes a set of internal 
state variables of the configurational subsystem. In the 
STZ context these are taken as concentrations of shear 
transformation zones in different orientations. For our 
application to SGR, we will take A as the distribution 
P[E, I). Note that because this distribution already con- 
tains information on all local elastic strains, we do not 
include separately in Uc a state variable for strain as was 
done in [13,111. 

The thermodynamic analysis is based on the first and 
second laws. The first law, i.e. energy conservation, takes 
the following form: 



Uc + Uk + Ur = ^(77 



(3) 



The r.h.s. is the external work rate that arises from 
shearing a system of volume V at shear rate 7 against 
a shear stress a. In terms of the subsystem entropies 
Sc, Sk and Sr, and the corresponding temperatures 
X — dUc /dSc\f^ and 9 = dUK/dSx, the first law can be 
rewritten as 



xSc + 0Sk 



Ur - Vaj + A ^^'^ 



OA 



= 



(4) 



Sc 



We only write one A-derivative term here, as a shorthand 
for a sum over the relevant derivatives w.r.t. the state 
variables gathered in A. 

The second law states that total entropy must increase, 
i.e. 

Sc + Sk + Sr>0 (5) 
One can eliminate Sc from this using the first law and 
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write Sr = Ur/Or to get 



1 - ^ I C^i? - ( 1 - 



esK 



-Ur 



where 



W = Vaj - A 



dUc 



dA 



W>0 
(6) 



(7) 



Sc 



Because time variations in Ur, Sk and A are in principle 
arbitrary, it is plausible to argue [l7[ that the three terms 
in ([6]) should be separately non-negative. To allow easier 
comparison with [16- 18], we follow this reasoning here. 
However, it turns out that this separate non- negativity 
is sufficient but not necessary for ([6]) : in Appendix |^ we 
present a less restrictive but still sufficient set of condi- 
tions for non-negative entropy production. Setting that 
aside and assuming separate non-negativity of all three 
terms, one obtains the conditions W > Q and 



-B( 



es 



K 



■^Ur 



'A[e-x) (8) 



where A and B are positive but can depend on the state 
of the system, e.g. via the various temperatures. Our no- 
tation for A follows Ref. [ll|, and we define B by analogy; 
the coefficients in [T6l - [l8| differ by factors of temperature. 
To simplify further, one assumes strong coupling between 
the thermal reservoir and the fast subsystem. In the cor- 
responding limit of large B, one then obtains 9 = 9r, but 
this equality is approached in such a way that the prod- 
uct B{9r — 9) stays finite. In fact it is only in this limit 
that ([U leads to dS]), as discussed further in Appendix El 
One can now write the first law as 



xSc^~{OSk + Ur) + W 



(9) 



Because of the assumption of strong coupling between the 
reservoir and the fast subsystem, one can exploit 9 = 9r 
to write the first term on the r.h.s. as —{9Sk + j-Ur) — 
A{9-x).so 



xSc^W + A{9 - x) 



(10) 



The upshot of this analysis [l7|, |l8| is thus twofold. 
Firstly, one obtains an equation of motion for the config- 
urational entropy; this also then determines the dynam- 
ics of the non-equilibrium temperature x as discussed in 
Sec. I VII The second result is the constraint W > Q. In 
([TU]). W is the dissipation rate, i.e. the difference between 
the external work (rate) and the (rate of) increase in the 
free energy of the A-degrees of freedom. 

To go beyond the formal expression for W above re- 
quires further assumptions on the relation between the 
state variables A and the other configurational degrees 
of freedom [T6l - [l8| . The slow subsystem as a whole is de- 
scribed via Uc{Sc, A) = Ua{A) + Ui{Sc - Sk{A)) where 
Ui{Si) captures all the slow degrees of freedom other 
than A. Then 



c 



dA 



Sc 



OUa _ dUi OSa 
~dA ~ dSi OA 



(11) 



One can now simplify using x — dUc/dSc\A — 
dUi/dSi. Assuming further that the number of inter- 
nal state variables in A is much smaller than the overall 
number of slow degrees of freedom fl^ , one has Sa ^ Sc 
and X can be regarded as a function of this vast majority 
of configurational degrees only. The final expression for 
the dissipation rate is then 



OSa 

- dA 



(12) 



and this shows that external work that is not dissipated is 
indeed stored in the free energy Ua — X^a of the internal 
state variables A. 



IV. APPLICATION TO THE SGR MODEL 

As mentioned earlier, in applying the above formalism 
to SGR our approach we take the distribution P{E, I) 
as the set of internal state variables A. The arguments 
above rely on the number of degrees of freedom in A 
being much smaller than the overall number of degrees 
of freedom in the configurational subsystem. This can 
be ensured by initially considering, instead of P{E,l), 
the total probability this distribution assigns to each of a 
finite number of bins covering the (E, I) plane. The num- 
ber of bins is then sent to infinity more slowly than the 
system volume V in the thermodynamic limit of large V. 
We assume below that this limit has been taken appro- 
priately and write directly the version of the calculations 
in the large system limit, i.e. without binning. 

Clearly for SGR one should choose as the internal en- 
ergy associated with A = P{E, I), 



UA{A)^Ne / dEdlP{E,l)u{E,l) 



(13) 



where N,. — V/ve is the number of SGR elements with 
volume Ve each, and u{E,l) is some function of the lo- 
cal yield energy and strain. The corresponding choice of 
entropy is less obvious. A natural proposal is S'a(A) = 
-TVe / dE dl P{E, l)[\iiP{E, 0-1]. However, it turns out 
that this choice does not give a thermodynamically con- 
sistent interpretation of the SGR model, essentially be- 
cause every yield event resets I — and so loses an infinite 
amount of entropy. We argue instead that, because the 
Z-dynamics is deterministic apart from jumps when rear- 
rangements occur, these degrees of freedom are "slaved" . 
Hence we should only consider the entropy of the distri- 
bution across yield energies. We also allow for a prior 
distribution (or density of states) R{E) in the entropy, 
and take: 



S'a(A) = -Ne / dEdlP{E,l) 



\ P{E) 
In— ^4 - 1 



(14) 



where P{E) = J dlP{EJ). The argument that the en- 
tropy should be independent of the /-degrees of freedom is 
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somewhat analogous to the observation made in STZ that 
in order to derive the standard STZ equations of motion 
one has to assume that the shear transformation zone 
orientation makes no contribution to the entropy [llj. 

To ensure that SGR fits into the thermodynamic 
framework described above, we need to show that with 
the above choices for J7a and Sa, one has V7 > al- 
ways, which means that the rate at which external work 
is performed is never less than the rate at which free en- 
ergy is stored reversibly in the internal state variables 
A = P{E,l). 

The inequality W > has to hold at all times t, so we 
drop all time arguments in the following. If we unpack 
our shorthand notation, the contribution to W in (fT2|) 
from the change of P{E, I) becomes 



A 



d{UA - xSa) 



dEdlPiE,l)^^-g^ (15) 



OA ~ J ' ' ' 5P{E,l) 

where the functional derivative on the right hand side is 



SP{EJ) 



u{E, I) + xlri 



P(E) 



R{E) 



(16) 



In combination with the external work, this gives for the 
dissipation rate 



P{E) 



W ^Vai-N,J dEdlP{E,l) , ^-^^^^ 

(17) 

We now want to check that W > 0, for the SGR equation 
of motion ([T]). Inserting the latter into the expression 
above and also expressing cr as an average over P{E, I) 
from (l2l) shows 



— = I dEdl {'^kvelP + idiP 



Xln 



P{E) 
R{E) 



{TP - Yp) 



Xln 



P{E) 
R{E) 



(18) 



where P, p, u are short for P{E,l), p{E,l), u{E,l) re- 
spectively. Now integrate by parts over I: 

^ = j dEdli^^PikvJ-diu) 



X{TP-Yp)\n- 



R{E) i 



(19) 



Since the sign of 7 is arbitrary, and so is the shape of 
P{E,l), the first term must vanish pointwise, i.e. diu = 
kvel- This of course makes sense: the stress contribution 
of each element is just the strain derivative of its internal 
energy density. So u{E, I) — kveP/2 is a function of E 
only, and from the meaning of the yield energy we expect 
this function to be —E itself so that u{E, I) — kvj? /2 — 
E. 

With the flow term thereby eliminated. 



-^^fdEdl fl^-pVn^^^ 
NeXY J \Y ' J R{E) 



(20) 



We now specialize to the conventional SGR model where 
p[E,l) = p{E)S{l), T{E,l) = roe"/x, and the obvious 
prior to use in the definition of S\ is R{E) = p{E). Also, 
the quantity 



TT{E,l)=T{E,l)P{E,l)/Y 



(21) 



is a normalized distribution by definition of y . So wc can 
express W as 

dEp{E)\n ^^^1 > (22) 

To show that this is non-negative we now use a cross- 
entropy inequality of the form 



dzP{z)lii[Piz)/Q{z)] > - In / dzQ{z) (23) 



This is valid if P is a probability distribution and Q is a 
measure (i.e. a probability distribution but not necessar- 
ily normalized to 1). The second term in W is already of 
the form ([23|) . In the first term, we can get tt{E, I) inside 
the log by writing P[E)P{l\E) = P{E, I); here P{l\E) is 
the conditional distribution of / given E. This yields 



N.xy-I'^'^''^^'^^'''T.pimpiE) 

+ JdEp{E)ln-P^^^ 



E/xp(E) ^^^^ 
> -ln(ro/r) - In J dEe-^/^P{E) (25) 

f dEdle^''^-^'/^'^'^/>^P(E,l) , , 
= '- ^ JdEdle-^fxPiEJ) 

This shows that the SGR model is thermodynamically 
consistent in the sense that its equation of motion en- 
sures that the dissipation rate W is non-negative as it 
should be. Note that this argument worked only because 
in writing down the yield rate T{E, I) above we had al- 
ready identified the SGR effective temperature x with 
the thermodynamic temperature x of the slow degrees of 
freedom. This is an important conclusion, which we will 
develop in more detail in the next Section. 

We remark finally that the three inequalities used 
above to prove W > (two cross-entropy ones, and the 
last one that the average of exp[fcwe^^/(2x)] over any dis- 
tribution is > 1) all become equalities if (and only if) 
P{E,l) oc p{E)e^/^5{l). As is physically sensible, this is 
the (Boltzmann) steady state of the SGR model without 
flow, where W = Q. In all other situations, W will be 
positive. 



V. CONSTRAINTS ON MORE GENERAL 
SGR-LIKE MODELS 

The natural question arising from the results of the 
previous Section is whether more general versions of the 
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SGR equation of motion are thermodynamically consis- 
tent in the same sense. 

A. Deviations from local linear elasticity 

Consider first deviations from local linear elasticity, as 
contemplated e.g. in the tensorial version of SGR 21]. 
These are simple to incorporate: one could have u{E, I) — 
u{l) — E, with the stress then the average of u'{l)/ve over 
all elements. As long as u(l) has its global minimum 
u{0) = at Z = 0, the argument above goes through, and 
the model is thermodynamically consistent. 

B. More general yield rates r{E,l) 

One can also generalize the SGR model by allowing a 
more general dependence of the yield rates T{E, I) on the 
state of the local element, i.e. E and /. (This was done, 
for instance, in [iS] in the context of shear thickening 
models; we discuss the implications for that work in Sec- 
tion |Vll) To understand the behaviour of W in this case, 
we start from the still general expression (PO)) . choosing 
again R{E) = p{E) as the prior in the entropy: 



W 



NeXY 



dE dl Tr{E, I) In 



/ 



dE p{E) In ■ 



P{E) 
P{E) 



-E/x 



P{E) 



(27) 



As before, one can now write P{E) = P{E,l) / P{l\E) = 
Y'K{E,l)/\r{E,l)P{l\E)] so that the first term becomes 



j dE dln{E,l)\n 

= J dEdlTT{E,l)\n 

Similarly, the second term is 
PiE) 



P{E) 



PiE) 

YTT{EJ)e''/^ 
T{E,l)P[l\E)p{E) 



(28) 



dEp{E) In 



e-E/xp{E) 



^p^/pr/ls^wpn nE,l)P{l\E)piE) 
dEdlP{l\E)p{E)\n ^-E/xYn{E,l) ^ ^ 

Putting these two together and isolating the cross- 
entropy terms gives 



W 



NexY 



= j dEdl [n{E,l)-P{l\E)p{E)]\n 



dEdlTT{E,l)\n 



dEdlP{l\E)p{E) In 



e"/x 

WT) 

T{E,l) 

e-Elx 



<EJ) 
P{l\E)p{E) 



(30) 



To say something about necessary conditions on r(£^,0) 
first, consider P{l\E) = 5{l), so that tt{E,1) oc 



r{E,0)P{E)6{l). Choose a P{E) such that this is 
p{E)[l+r]r{E)]S{l) with some small t] and a so far unspec- 
ified function r{E) satisfying / dE p{E)r{E) — 0. Then 
the first contribution to W is O^rf'). The second and 
third combine to 



■q / dE p{E)r{E)\n- 



~E/x 



Unless the argument of the logarithm is a constant, one 
can arrange the sign of r{E) so as to make the integral 
positive, and thus W negative to leading order in 77. To 
avoid this contradiction to > 0, r{E, 0) must be pro- 
portional to e-^/>^, say T{E,0) = roe"-^/^. This ar- 
gument shows that for thermodynamic consistency, the 
effective temperature x of SGR must be equal to the non- 
equilibrium thermodynamic temperature x of the config- 
urational subsystem, as promised in the previous Section. 

Turning next to ^-dependence of r(i?,Z), we can write 
this as T{E,l) = roe-^/^g{E,l) with giE,0) = 1, and 
hence 

W /",„.._„ T^iEJ) 



NeXY 



dEdl [tt{E,1) ~ P{l\E)p{E)]\n 



dE dlTT{E,l)\n- 



,fet-e/V(2x) 

9{E,l) 



dE dlP{l\E)p{E) \ng{E,l) 



P{l\E)p{E) 



(32) 



The first term here is of cross-entropy form, while the 
second and third term are (pointwise) non-negative if 1 < 
g{E,l) < e^"""^^ I'y'^^) . In summary, a sufficient condition 
for > is 

r(£;, = Vne-^l^g{E, I) with 1 < g{E, I) < e'="'=''/(2x) 

(33) 

The first part is also necessary as shown above. Whether 
the second part is likewise necessary is not clear, i.e. there 
may be other g{E,l) that also keep > 0. But physi- 
cally the sufficient condition given here is broad enough: 
the yield rate should increase with P, i.e. it should never 
go below its Z = value for given E; on the other hand 
the yield rate should not go up more quickly with P than 
prescribed by the reduction in activation energy. If any- 
thing, one might expect the rate to level off when the 
barrier —u = E — kvel'^ becomes small or even negative, 
for example 



or 



r(£;,/) =romin(e^/'^,l) 



r(ii;,o = ro/(i + e-"/'^) 



(34) 



(35) 



It is clear from the above arguments that one could also 
allow a shear-rate dependence of the attempt rate Fq, e.g. 
replacing Fo by (Fg -I- 7^)^^^ in the spirit of recent pro- 
posals [22| on how the relaxation time varies with shear 
rate. These more general forms can of course be com- 
bined with deviations from local linear elasticity as in 
the previous subsection, with the upper limit on g{E, I) 
then e"(')/'^ instead of e'=''=''/(2x). 
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C. Frustration 

We had previously considered [l3] the possibihty that 
the local strain I might not relax fully to zero after a yield 
event because disordered stresses from the environment 
prevent this, and termed this effect 'frustration'. Math- 
ematically, it would correspond to p{E,l) — p{E)r{l\E) 
with r{l\E) ^ 5{l). For an SGR model extended in this 
way it is doubtful that an argument can still be made 
that I should be left out of the entropy S'a as / no longer 
has deterministic dynamics slaved to the last yield time 
of each element. If one nevertheless leaves S'a as before, 
then one can show that the extended model does not 
guarantee > 0, i.e. there are choices for P{E, I) which 
lead to < 0. It is possible that the model would still 
be thermodynamically consistent because such distribu- 
tions might not be accessible by the dynamics starting 
from reasonable initial conditions, but we have not been 
able to establish this. 



VI. SUMMARY AND DISCUSSION 

With plausible expressions for the entropy and inter- 
nal energy of the P{E, Z)-degrees of freedom, we have 
been able to show that the SGR model admits an in- 
terpretation consistent with the non-equilibrium thermo- 
dynamics framework of Bouchbinder and Langer. This 
interpretation requires that the effective temperature x, 
postulated phenomenologically within the SGR model, 
is in fact the thermodynamic temperature x of a con- 
figurational subsystem describing mesoscopic degrees of 
freedom. (These degrees of freedom have slow energy 
transfer to the thermal bath of 'fast modes' and hence 
their temperature is a nonequilibrium one.) 

We believe this represents significant progress in under- 
standing the conceptual foundation of the SGR model. 
The original empirical motivation for introducing x 
was (and remains) intuitively reasonable: an effective 
noise temperature arising ultimately from yield events 
throughout the sample. However a more solid formal 
grounding is now provided by the link to a thermody- 
namically defined temperature that governs the explo- 
ration of the system's inherent structures, through the 
dynamics of its mesoscopic degrees of freedom. These 
degrees of freedom are taken to be in thermal equilib- 
rium with each other, but not with the rest of the sys- 
tem. If one accepts that, the thermodynamic interpre- 
tation gives new insights by showing, for example, that 
X should be viewed as intrinsic to the state of the sys- 
tem, and thereby to its history, as opposed to being set 
directly by the current rate of external driving. (The 
latter would make x-mediated effects second order, and 
hence negligible, within linear response theory [23j.) 

Of course, the consistency of the thermodynamic in- 
terpretation of X does not make this the only possible 
interpretation of the SGR model. It would be perfectly 
defensible to believe that, in many soft glasses, meso- 



scopic degrees of freedom do not equilibrate sufficiently 
even among themselves for their non-equilibrium temper- 
ature to ever be defined. In that case, however, it re- 
mains useful to know that the empirical choice made for 
the noise-driven dynamics in the standard formulation 
of SGR is exactly equivalent to having such a temper- 
ature. Anyone convinced that such a temperature does 
not exist would probably wish to add terms to break this 
equivalence - which amounts to an implicit assumption 
of detailed balance among slow modes in the quiescent 
state. 

We have further shown that the thermodynamic con- 
sistency extends from the original SGR model to more 
general versions that allow for nonlinear local elasticity. 
Likewise, a fairly broad class of yield rates r(£', I) is al- 
lowed. This includes most physically plausible choices, 
but always maintaining the thermally (x-)activated de- 
pendence on yield energy. 

Beyond checking consistency, the thermodynamic 
framework gives an equation of motion (jlOp for the en- 
tropy Sc of the configurational subsystem; see the Ap- 
pendix for a generalization of this. As discussed in 
one can convert Eq. (|10l) to an equation of motion for 
the effective temperature (a more involved scenario is 
discussed in ^^). To see this, recall the decomposi- 
tion we assumed for the energy of the slow subsystem, 
Uc{ScA) = Uk{K) + Ui{Sc - Sa{K)). The effective 
temperature x — dUc /dSc\/^ — dUi/dSi depends then 
on Si = Sc — Sa- Assuming as before that S'a <C Sc 
because A only contains a vanishing fraction of the slow 
degrees of freedom, x becomes a function of Si « Sc 
only. Conversely, x then fully determines Sc, so that 
Sc = {dSc/dx)x and Eq. ^ becomes 



Cfx ^W + A{e-x) 



(36) 



Here C^(x) is an effective heat capacity at constant 
volume, and A{x, 0) is as before positive but otherwise 
unspecified; in principle, A could also depend on other 
quantities like 7 and possibly even P{E,l). What is in- 
teresting, however, is that the form of the driving term 
W is not negotiable; this is from (|20)) 



W 



J dE dlT{E,l)P{E,l)\n 



P{E) 



- Y / dEp{E)\n 



P{E) 
P{E) 



(37) 



While not obvious from this expression, it is clear from 
the definition Q that in steady state, where A = P{E, I) 
does not change, W = Va'^ is just the external work rate. 
This is the form that was assumed recently in exploring 
the extension of the SGR approach to understand shear 
banding: it was called 'Model 1' in [l5|. Steady state 
flow curves of shear stress tr versus 7 as predicted from 
p6| would therefore look exactly like those found in [l^ . 
On the other hand, 'Model 2' from that paper, where 
the driving term for x is taken as proportional to the 
yield rate Y ^ is not thermodynamically consistent. We 
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note that even with model 1, Eq. (|37|) predicts different 
transient behaviour, because the driving term W gener- 
ally differs from the external work rate Vcr-y. It would 
be interesting to explore in future work what effects this 
thermodynamically motivated form for W has on the pre- 
dicted shear banding dynamics. 

We can similarly assess the thermodynamic consis- 
tency of two extensions of the SGR model that were 
proposed some time ago to model shear thickening ef- 
fects [l^. That paper considered two models. In one, 
the effective temperature a; = x is taken as a function 
x{a) of the overall stress. The second model, in a more 
radical departure from the standard SGR picture, takes 
x{l) as a function of the local strain in each element. The 
first case can feasibly be accommodated within the ther- 
modynamic framework, at least if we consider situations 
where the effective temperature x, in its time evolution 
according to ([36|). has reached a steady state. Assuming 
that 9 <ti X for simplicity, so that the effective tempera- 
ture is much larger than the bath temperature, the steady 
state condition is x = W/A = Vaj/A. Given that A is 
only restricted to be positive but allowed to depend on 
e.g. 7 and a, there are certainly choices for this depen- 
dence that would give a steady state x depending (only) 
on a. 

The second model in [l^, where a; = a::(^), is different. 
Having a temperature that differs from one element to 
the next is a significant departure from the idea of a 
single effective temperature for the slow configurational 
subsystem. Consistent with this, one can check that a 
generically chosen x{l) will violate the conditions (|33| . 
This does not of itself prove W < 0, but we have checked 
that there are indeed conditions under which negative 
W does arise, e.g. for the specific (stepwise decreasing) 
form of x(l) considered in ^19;]. Accordingly this second 
model of [l9| is not consistent with a thermodynamic 
interpretation of the effective temperature as described 
here. 

Finally, we should note again that our thermodynamic 
viewpoint for SGR was developed above by closely fol- 
lowing the route first taken for STZ by Bouchbinder and 
Langer [T6l-[T8j . These authors used the second law to in- 
fer conditions on the dynamics of macroscopic variables. 
However these conditions, while both plausible and suf- 
ficient, are not in fact necessary: in the Appendix we de- 
scribe a less restrictive, but still sufhcient, set. The true 
value of thermodynamic consistency as a constraint on 
mesoscopic rheological models will not be known with- 
out a set of necessary, as well as sufficient, conditions 
for compliance of such models with the second law. It 
is possible that an explicitly model-independent formal- 
ism, such as that developed in [2^1, could prove useful in 
identifying such conditions. 

A further interesting direction for future work will be 
to explore thermodynamic effects on the aging behaviour 
predicted by the SGR model. An important feature of 
the model is that aging effects are present [H, even 
when X is constant, as long as its value lies below the 



glass transition a,t x = Eq. It is clear that allowing x 
to evolve dynamically according to (|36p in line with the 
thermodynamic picture will contribute additional aging 
effects. However, what quantitative form these take can 
only be fully assessed once physically motivated param- 
eter dependences have been established for the effective 
heat capacity and the coefficient A determining the 
rate of heat transfer between the slow and fast subsys- 
tems. 
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Appendix A: More general conditions for second law 

We have in the main part of the paper adopted the 
argument of that the three terms in © must be sep- 
arately non-negative. This is a sufficient condition for 
the overall sum to be non-negative and is also physically 
plausible. However it is not mathematically necessary. 
We discuss here a somewhat more general set of condi- 
tions that are also sufficient. These admit the thermody- 
namic consistency of a wider class of models than allowed 
by the conditions of [13] as used above. 

If we write the first law as in ^ and insert Ur = GrSu, 
we obtain 

xSc + OSk + OrSr = W (Al) 
This needs to be satisfied together with the second law 



Sc + Sk + Sr > 



(A2) 



With the reversible degrees of freedom (A) of the slow 
system already treated separately in W, these two equa- 
tions are completely symmetric in the roles played by the 
slow, fast and reservoir subsystems. A more general set 
of conditions for the second law to hold then suggests 
itself as ly > and 

xSc = acW + A{9 - x) + C{0r - x) (A3) 
OSk = axW + A{x -9)+ 3(6^ - 6) (A4) 
OrSr = urW + B(e - du) + C(x - Or) (A5) 

where A, B, C, ac, otK-, cur are all non-negative and 
ac + ctK+OiR — 1- The last condition is enough to ensure 
that the first law is satisfied, as the heat flow terms all 
cancel in pairs. For the total rate of entropy change one 
gets 



Sk + Si, 



W 



ac OLK_ ail 

X 0R 



A 



{x-0? 

X9 



^iO^On)\^ix-ORr 
xOr 
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which is evidently non-negative. The more restrictive 
condition of [iTj corresponds essentially to = ct^ — 
C = 0. One has then Un = OrSr = -B{eR - 9) as in 
(©, and 

ffl OR 

This generalizes the second equality of (jS]) to the case 
where Q ^ Br. In the limit where B becomes large at 
constant B^Qb. ~ 0), the term on the r.h.s. proportional 
to B goes to zero so ([8]) is retrieved. That a term pro- 
portional to B{9r — 9)^ should be present in general here 
is clear because 9Sk + [9/9r)Ur = 9{Sk + Sr) and this 
rate of entropy change should have a non-zero contribu- 
tion when there is heat transfer at a finite rate between 
the reservoir and the kinetic-vibrational subsystem. 

Returning to our more general conditions above, the 
effect of the coefficient C is probably unimportant so long 
as we are indeed interested only in the limit 9 — 9r 
because of strong thermal coupling to the reservoir: one 
can then just combine A and C. 

The effect of the a coefficients is more subtle. In soft 
materials, where the reservoir should include the solvent 
degrees of freedom, it seems plausible that a substantial 
amount of the work performed on the system (minus the 
amount stored reversibly in A, which gives W) could be 
dissipated directly in the solvent. For instance, in a dense 
emulsion, work could be dissipated directly in the solvent 
films between droplets without ever passing through the 



A degrees of freedom. In soft glasses, therefore, aR could 
be substantial, and correspondingly ac < 1- An alterna- 
tive strategy might be to exclude this contribution from 
the work input, by separating off a global Newtonian sol- 
vent contribution. However, as the example of thin films 
between emulsion droplets shows, the mesoscopic and sol- 
vent contributions cannot be treated as simply additive: 
in principle the droplet organization controls both the 
thicknesses of films and the shear rates within them. It 
therefore seems advisable not to assume that kind of sep- 
aration within a thermodynamic description, at least in 
the case of soft glasses. 

Despite the above generalizations the condition that 
> remains in place. So far we have not been able to 
show mathematically that W > is necessary (or indeed 
that there are not broader conditions on the heat flows 
that would satisfy the second law), although physically 
it is certainly plausible that the dissipation rate should 
be non-negative. The main possible generalization aris- 
ing from the arguments above is then in the equation of 
motion for which would become (replacing ([36])) 

Cfx = acW + Ai9 - x) (A8) 

Here one now has an additional undetermined factor 
< ac < 1 in front of the driving term W. One 
could speculate whether there is a difference between e.g. 
metallic and soft glasses in the size of ac, which could 
be of order unity for the former and small for the latter. 



For a recent overview, see e.g. the Topical Issue on the [141 S. M. Fielding , P. Sollich, and M. E. Gates, J. Rheol 



Physics of Glasses, Eur. Phys. J. E, 43(9), 2011. 



J. M. Brader, T. VoiRtmann, M. E. Gates, and M. Fuchs 



Phys. Rev. Lett., 98, 058301 (2007)[ 

M. Brader, M. E. Gates, and M. Fuchs, j'hys. Rev 



44, 323 (2000)1 

1 151 S. M. Fielding: M. E. Gates, and P. Sollich, ^oft Matter 



Lett., 101, 138301 (2008) 



5, 2378 (2009) . 
1 16 1 E. Bouchbinder and J. S. 



J. M. Brader. T. Voigtmann. M. Fuchs. R. G. Larson. 



031131 (2009) . 
[17j E. Bouchbind er and J. 



(2009H 



and M. E. Gates, |PNAS, 106 15186 

P. Hebraud and F. Legueux. |Phvs. Rev. Lett.. 81. 2934 



(1998) 



031132 (2009) 



[18] E. Bouchbind er and J. S 



Rev. E, 66, 051501 (2002) 



031133 (2009) 



Phys. 


Rev. 


E, 


80, 




Phys. 


Rev. 


E, 


80, 




Phys. 


Rev. 


E, 


80, 



G. Picard. A. Aidari. L. Bocguet. and F. Legueux. Phys. [19] D. A. Head. A. Aj dari. and M. E. Gates, Phys. Rev. E 



L. Bocguet. A. Goli n. and A. Ajdari, Phys. Rev. Lett 



103, 036001 (2009) 



64, 061509 (2001) 



O. Pouligucn and Y. Forterre, 



A , 367, 5091 



M. Talamali, V. Petaja, D. Vande mbroucg , and S. Roux 



Phj^s. Rev. E, 84, 016115 (2011 



'2009J 



Phil. Trans. R. Soc. Lond 



[20] M. Tsamados, A. Barra, a nd P. Sollich. unpublished 
[21] M. E. Gates and P. Sollich, |j. Rheol., 48, 193 (2004") . 
[22] H. M. Wyss, K. Miyazaki, J . Mattsson, Z. B. Hu. D. R, 



Reichm an, and D. A. Weitz, Phys. Rev. Lett., 98, 238303 



(2007) 



V. Mansard, A. Golin, P. G hauduri, and L. Bocquet 
Soft Matter, 7, 5524 (25TT| 



1 23 1 E. Bouchbinde r and J. S. Langer, Phys. Rev. Lett., 106 



M. L. Falk and J. S. Lan ger, Ann. Rev. Gond. Matt 



Phys., Vol 2, 2, 353 (2011 



148301 (2011)1 

|24| E. Bouchbinder and J. S. Langer, ^oft Matter, 6, 3065 



(2010) 



Phys. Rev. Lett.. 78, 2020 fl997| 


P. Sollich, 


Phys. Rev. E, 58, 738 (1998) 



[25] H. G. Ottinger, Beyond equilibrium thermodynamics 
(Wiley, New Jersey, 2005). 



